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^ \ Abstract. The BKP hierarchy has a two- component analogue (the 2-BKP hierarchy). Dis- 

^ ■ persionlcss limit of this multi-component hierarchy is considered on the level of the r-func- 

tion. The so called dispersionless Hirota equations arc obtained from the Hirota equations 
of the r-function. These dispersionless Hirota equations turn out to be equivalent to a sys- 
tem of Hamilton-Jacobi equations. Other relevant equations, in particular, dispersionless 
Lax equations, can be derived from these fundamental equations. For comparison, another 
C/j , approach based on auxiliary linear equations is also presented. 
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q ! 1 Introduction 

O, 

Among many variants of the KP hierarchy PJ |2j, the BKP hierarchy [HI EJ El IE] enjoys a dis- 

\q , tinguished status. First of all, both the KP and BKP hierarchies are formulated by a single 

r-function [2j. The r-function of both hierarchies has a fermionic representation in terms of 

a one-component free fermion system [I]. The only difference lies in the type of fermions 

(charged or neutral). In the context of algebraic geometry, the KP hierarchy characterizes Ja- 

cobi varieties, and the BKP hierarchy corresponds to Prym varieties [51[71IHj. As regards the 

Lax formalism, both hierarchies are formulated by scalar pseudo-differential operators j^j. From 

S^ . this point of view, the BKP hierarchy is a reduction of the KP hierarchy. This relation carries 

5-j | over to the dispersionless limit [H]- 

The subject of this paper stems from the last topic, namely, the dispersionless BKP hierar- 
chy. As for the general background on dispersionless limit of integrable hierarchies, we refer to 
Takasaki and Takebe's review |l(Jj . 

Recently, Bogdanov and Konopelchenko studied the dispersionless BKP hierarchy from a new 
approach [TT] . Their approach is based on the 9-dressing method |12l 1131 1141 I15j , which had 
been successfully applied to other dispersionless integrable hierarchies including the (genus zero) 
universal Whitham hierarchy ^SJ- Employing this method, they obtained an analogue of the so 
called "dispersionless Hirota equations" for the BKP hierarchy, which turns out to be slightly 
different from those of the KP and Toda hierarchies [El H3 HH1 HH1 E3 EU ■ 

The aim of this paper is to re-examine this issue along the lines of approach originally taken 
for the dispersionless Hirota equations |1UI I17j . Namely, we now start from the Hirota equations 
for the r-function itself rather than its dispersionless limit. Moreover, we consider the two- 
component BKP (2-BKP) hierarchy that includes two copies of the usual one-component BKP 
(1-BKP) hierarchy as a subset. 
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Unlike the 1-BKP hierarchy, the 2-BKP hierarchy is not very well known in the literature (see 
Kac and van de Leur's paper (22] for a comprehensive account of this kind of multi-component 
integrable hierarchies). The 2-BKP hierarchy was first discovered as a system of Hirota equa- 
tions |22j- A Lax representation was later proposed in an algebro-geometric framework 0IH], 
which revealed that this system is closely related to two-dimensional finite-gap Schrodinger 
operators and their isospectral deformations (the Novikov-Veselov hierarchy) |24j . 

The 2-BKP hierarchy is rather special among multi-component integrable hierarchies, because 
it has a scalar Lax formalism as Shiota [Jj and Krichever [B] presented. In this respect, the 2- 
BKP hierarchy resembles the Toda hierarchy. The Toda hierarchy can be interpreted as the two- 
component KP (2-KP) hierarchy, though it has a scalar Lax formalism (in terms of difference 
operators) [22]. The scalar Lax formalism of the 2-BKP hierarchy, however, has a new feature. 
Namely, it has two spatial dimensions unlike the one-dimensional setting for the KP, Toda and 
1-BKP hierarchies. This is the origin of all complexities that we shall encounter. 

Our main concern lies in the structure of dispersionless Hirota equations and some other 
relevant equations in such a multi-component setting. We first derive a set of "differential Fay 
identities" from the Hirota equations (more precisely, their generating functional form) . We then 
apply the standard prescription of dispersionless limit |10(ll7j (which is a kind of "quasi-classical 
limit") to these equations. The outcome is a set of equations that may be called dispersionless 
Hirota equations. Whereas part of them coincide with the equations presented by Bogdanov 
and Konopelchenko ^J, other equations are seemingly new and specific to the two-component 
case. Moreover, we can correctly derive other relevant equations of the dispersionless hierarchy 
such as Hamilton-Jacobi equations, Lax equations, etc. We shall see that two distinct Poisson 
brackets show up in the formulation of dispersionless Lax equations. They are manifestation of 
presence of two spatial directions mentioned above. 

Let us mention here that Konopelchenko and Moro have studied the dispersionless limit of 
the Novikov-Veselov hierarchy by means of (9-dressing |26j . Our results should be, in principle, 
related to theirs. 

This paper is organized as follows. Sections 2 and 3 are devoted to a brief account of the tau 
function and equations satisfied by the tau functions. In Sections 4 and 5, we derive the diffe- 
rential Fay identities and the dispersionless Hirota equations. In Sections 6 and 7, we examine 
their implications to the Hamilton-Jacobi equations and the dispersionless Lax equation. In 
Section 8, these results are compared with the conventional approach based on auxiliary linear 
equations. 

2 Tau functions and bilinear equations 

The 2-BKP hierarchy has two sets of time variables £i, £3, . . . , t2n+ii • • • and i 1 , i 3 , ■ ■ ■ , t2n+i: ■ ■ ■ 
with odd indices 1 . For notational convenience, We collect these variables into two oo-dimensional 
vectors 

t = (h,t 3 , . . .), t = (*i, *3, . . .). 

The r-functions r = r(t,t) of this hierarchy are characterized by the fundamental bilinear 
equation [201221 



/ ^-/^'-'Mt' - 21*- 1 ], i')r(t + 2[z-\ t) 



Z e^'-M) T (t', t' - 2{z- 1 ])r(t, t + 2[z' 1 ]) (2.1) 



2-kiz 



throughout this paper, the bar "~" is used to create a new symbol, meaning nothing like complex conjugation. 
For instance, t n and £„ are understood to be independent variables. 



Dispersionless Hirota Equations of Two-Component BKP Hierarchy 3 

that holds for arbitrary values of i! , t, i and t. Both sides of this equation are understood to 
be contour integrals along the circle \z\ = R with sufficiently large radius R and anti-clockwise 
orientation, i.e., 

Alternatively, these integrals may be thought of as a purely algebraic operation 
for (formal or convergent) Laurent series. £(£, z) and [z~ 1 ] denote the BKP version 

oo i _3 _ 2n— 1 \ 

e(M = £w^ + \ [z- 1 ] = <z-\ z —,...,^- TT ,...) 



n=0 



of the usual notation for the KP hierarchy. 

We can rewrite this bilinear equation to an infinite number of Hirota equations as follows. 
By substituting 

t n ► t n x ni t n > t n -\- x n , t n > t n x n , t n > t n -\- x ni 

(|2.1|) transforms to 

I —e 2 ^ x ^T(t -x- 2[z~\ i - x)r{t + x + 2[z~\ t + x) 

= I —e 2 ^ z \(t -x,t-x- 2[z- 1 ])T(t + x,t + x + 2[z' 1 ]) , 

where x = (x±,X2, . . .) and x = (xi,x~3, ■ ■ .). If we use the Hirota bilinear operators D t2n+1 
and Dj 1 that act on the product of two functions f(t,i), g(t,t) as 

A 2m+1 /(M) ■ g(t,t) = d t2n+1 f(t,t) ■ g(t,t) - f(t,t) ■ d t2n+1 g(t,t), 
A 2m+1 /(*,*) • <?(*,*) = df 2n+1 f(t,i) ■ g(t,t) - f(t,t) ■ dj 2n+1 g(t,i), 

the product of the r-functions in the integral can be expressed as 

r(t + x + 2[z~\ t + x)T(t-x- 2[z~\ t-x) 

= exp(£(2£> t , z) + x ■ D t + x ■ D$r{t, t) ■ r(t, t), 
r(t + x,i + x + 2[z~ 1 ])r(t -x,t-x- 2[z~ 1 ]) 

= exp(Z(2D- t , z)+t-D t + t- D- t )r(t, i) ■ r(t, t), 



where 



1 1 

D t = [ D tl , - A 3 , • • • , 2n + l Dt2n+1 ' 



1 1 

3 * 3 '-"' 2^+1 



D t - [ D h , q A 3 , • ■ • , Qm , -, A 2n+1 i • • • 



X D t = ^ X 2 n+lD t2n+1 , X-D t = Y^ X 2n +l A 2 „+i ■ 

n=0 n=0 
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Let us introduce the polynomials qj(t), j = 0, 1, 2, . . . defined by the generating function 

oo 
3=0 

We can thereby rewrite (|2.1[) as 

oo 

J2 9j(2x) qj (2D t ) exp{x D t + x- Di)r{t, t) ■ r(t,i) 

3=0 

oo 

= J^ qj {2x) qj (2bi) exp(:r £> t + i • £>t)r(t, t) • r(t, i). 

3=0 

Taylor expansion of this equation with respect to x and x gives an infinite number of Hirota 
equations. In other words, this is a generating functional form of those Hirota equations. 

The 2-BKP hierarchy contains two copies of the 1-BKP hierarchy in the t- and t-sectors. 
These subhierarchies show up by setting i = i or t! = t. If we set t =i, (|2.1[) reduces to 

f 2^ eC(t '~ Mr(t ' " 2[Z ~ 1] ' *) r( * + 2[Z ~ 1] ' *) = r( *'' * )r( *' * } - 
Similarly, setting t' = t yields 

I -^-/ (1 '~ 1,Z)t ^ l ' ~ 2[* -1 ]M*> * + 2[^ -1 ]) = r(t, t')r(t, t). 
These equations are exactly the bilinear equations of the 1-BKP hierarchy [31 E] . 

3 Fermionic representation of tau functions 

The r-functions of the 1-BKP and 2-BKP hierarchies have a fermionic representation in terms 
of neutral fermions |H 1221125] . Let ipj, j € Z, denote the generators of a Clifford algebra with 
the anti-commutation relations 

VPji¥k]+ = (-l) J 5 j+kfi , 

and |0) and (0| the vacuum states of the Fock and dual Fock spaces that are annihilated by half 
of these generators ( "annihilation operators" ) as 

(0\ifj = for j > 0, ifj\0) = for j < 0. 

This is the fundamental setting of the neutral fermion system. Vacuum states of this system 
are two-dimensional. Namely, there are another pair of states |1) and (1| that are annihilated 
by the same set of annihilation operators. These vacuum states are interchanged by the action 
of V2ifo as 

y/2ipo\0) = \l), V2<po\l) = \0), (0\V2<p = (1\, (l\y/2<p = <0|. 

Now suppose that g is an invertible element of the Clifford algebra that generates a "Bogo- 
lioubov transformation" of the form 



dfkg 1 = ^2 a i kLp i 
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on the linear space spanned by 0^'s. A typical (generic) case is the exponential 



9 = exp I - ^2 a jkfjfk 

\ j,k>0 J 



of a fermionic bilinear form. Such an operator determines a r-function of the 2-BKP hierarchy 
as the vacuum expectation value 

r(t,t) = (0|e^V" 5 ®|0), 
where H{t) and H(t) are linear combinations 

oo oo 

H(t) = 2 J t2n+lH2n+l, H(t) = £ ^n+lH^n-l 

n=0 n=0 

of the "Hamiltonians" 

H 2 n+l = -^2(-iy +1 fjf-j-2n-l- 

Note that this resembles the fermionic representation of the r-function of the Toda hierarchy 

H3I2H1 

The fundamental bilinear equation (|2.1I) is a consequence of the algebraic relation 
^{- l ) 3l Pj9 ® <P-j9 = X^ -1 ^ 5 "^' ® W-i 

of the foregoing operators and the "boson- fermion correspondence" 

V2(l\e H ^ V (z)\V) = X(t, z)(0\e H ^\V), 
V^iWltpiz-^e-^ll) = X{i,z)(W\e- B{l) \0) 

that holds for arbitrary states (W\, \V) in the Fock and dual Fock spaces. <fi(z) and X(t,z) 
denote the free fermion field 

<p(z) = ^2<pj z3 

jez 
and the corresponding vertex operators 

oo _2n_l 

X(t,z) = e ^e-™^\ D(t,z) = £ ^^ . 

71=0 

The operator D(t, z) plays a fundamental role in the formulation of dispersionless Hirota equa- 
tions. 

4 Differential Fay identities 

Employing the method developed for the KP hierarchy jlUl I17j . we can derive four types of 
"differential Fay identities" . They are derived from (|2.1|) by the following procedure. 

1) Differentiate (|2.1j) by t'i and set 

t' = i + 2[A~ 1 ] + 2[^~ 1 ], t' = t. 
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2) Differentiate (|2~Tj) by t[ and set 

t' = t, t' = t + 2[A _1 ]+2[/i -1 ]. 

3) Differentiate (|2.1[) by t'i and set 

t' = t + 2[A~ 1 ], t' = i + 2\p,' 1 ]. 

4) Differentiate (|2.1|) by i[, and set 

t' = i + 2[A~ 1 ], t' = t + 2[iT 1 ]. 

A and /U are arbitrary complex constants that sit on the far side of the integral contour, namely, 
|A| > R and |/i| > R. 

1) In this case, we obtain a bilinear equation of the form 



/ 



dz (z + \)(z + im) f 



{zT(t + 2[\- l ] + 2[fi- L ]-2[z- L ],i) 



2iriz (z — X)(z — /i) 
+ (d tl r)(t + 2IA" 1 ] + 2[fi- 1 } - 2[z~\ t))r(t + 2[ z -\ t) 

= § 2^(^)(* + 2[A_1] + 2[ ^ 1] ' * " ^M*' * + 2 ^ -1 ])' 

where (9i 1 r)(t, t) denotes the ^-derivative of r(t,t). Also the well known formula 



00 ~ y 2n+l 



2 z ^+i i + z 

> = log 

^ 2ra + 1 & 1 - z 

n=0 

has been used to rewrite the exponential factor on the left hand side as 

gtftf-M) = e 5(2[A- 1 ]+2[ M - 1 ], 2 ) = (z + A)(z + /z) 

(z-\)(z + fj,y 

Since the integrands on both sides of this bilinear equation are meromorphic functions on the 
far side \z\ > R of the integration contour, the integrals become the sums of their residues at 
poles. Whereas poles of the integrand on the left hand side are located at z = A,/i, oo, the 
integrand on the right hand side has the only pole at z = oo. We thus obtain a somewhat 
messy bilinear equation for special values of the r-function and its ti-derivative. After some 
algebra, this bilinear equation boils down to the following relatively simple form that resembles 
the differential Fay identity of the KP hierarchy |1UL 117) : 



\-ti\ ' " T(* + 2[/i-i], I) 

A + „ - 3„ log T " + 2 ' A ''' + 2 fr"'- f>) T ," + f ,'' t 2 f '• ft' B - (4.1) 
' r((,t) ,/ T(( + 2[A-i], *)t(( + 2[m-'], t) 

Since these calculations are related to the t-sector only, the same equation holds for the tau 
function of the 1-BKP hierarchy as well. 

2) This case is essentially the same as case 1) except that the roles of t and t are interchanged. 
The final form of the bilinear equation reads 

A + M/\ « , r(t,t + 2[A- 1 ]) 



A — ji 



A — u — Of ■ log — — = F — ^ 
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A , , ~ , r(t, j + 2JA- 1 ] + 2[^ 1 ]) \ r(t, t + 2JA- 1 ] + 2[ / u- 1 ])r(t, t) 
-^A + M-^log ^ J T (t,t + 2[A-i])r(t,t + 2[^-i])- ^ 

3) In this case, we have a bilinear equation of the form 

+ (d tl r)(t + 2[A" 1 ] - 2[z~ 1 ], t + 2[^ 1 ]))r(t + 2IZ" 1 ], t) 

= f^-^(d tl r)(t + 2[X' 1 ],t + 2[^])T(t,t + 2[z^]). 

Again by residue calculus, we obtain the following analogue of the differential Fay identity: 
r(t + 2[A- 1 ], t) 



A - d tl log 



r(t, t + 2^" 1 ]) 



A - d tl log r(t + 2[A-i],t + 2^])\ rft + 2[A-i],i + 2^])r(t,^ 

1 S r(t,t) / r(t + 2[A- 1 ],t)r(t + 2[/ i - 1 ],t) l ; 

Unlike the previous two equations, this equation has no counterpart in the 1-BKP hierarchy. 

4) This case is parallel to case 3). The roles of t and i are simply interchanged. The final 
form of the bilinear equation reads 

t(M + 2[aT 1 ]) 
"-^ log T(t + 2[A-i],t) 

r(t + 2[A~ 1 ], t + 2[ / u~ 1 ]) \ r(t + 2[A~ 1 ], t + 2[ / u~ 1 ])r(t, j) 
M " t_1 g rjtj) jr(t + 2[A-i],t)r(t + 2M,tT (44) 

This equation, too, has no counterpart in the 1-BKP hierarchy. 

We have thus derived the four types of equations, (|4.1|) - (|4.4|1 . as an analogue of the differential 
Fay identity of the KP hierarchy. The latter is known to be equivalent to the full KP hierarchy 
jlUl Appendix B] . It seems likely that the same is also true for the 2-BKP hierarchy. Namely, we 
conjecture that the system of equations (|4.1|) — (|4.4|) will be equivalent to the 2-BKP hierarchy. 

The subsequent consideration may be thought of as partial evidence that supports this con- 
jecture. We shall examine the dispersionless limit of these equations, which will eventually turn 
out to be equivalent to the dispersionless limit of the 2-BKP hierarchy. 

5 Dispersionless Hirota equations 

As in the case of other integrable hierarchies JU], dispersionless limit is realized as a kind of 
"quasi-classical limit". To consider this limit in the language of r-function, we allow the r- 
function r = t(H, t, i) to depend on the "Planck constant" H and assume that the rescaled 
r-function 

Th (t,i) = T(h,h- 1 t,h~ 1 t) 

behaves as 

log Th (t, I) = h~ 2 F(t, i) + 0(h- v ) (5.1) 

in the classical limit h — > 0. The goal of the following consideration is to derive a set of differential 
equations for the scaling function F = F(t,i) from the differential Fay identities Q4.ljl - H4.4j) . 
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Upon rescaling the time variables as above, the derivative operators and the bracket symbols 
in (|4.1|) — (|4.4|) are also rescaled as 

d t2n+1 -»• hd t2n+1 , d hn+1 -^hdt 2n+1 , [A" 1 ] -> ^A" 1 ], [/i _1 ] -> % -1 ]. 

Bearing this in mind, let us examine the limit of these equations as h — > 0. 
1) Let us first consider (|4.1|) . By rescaling, this equation transforms to 

1 A - // - ftd^ log 



r ft (t + 2^[A- 1 ]+2%- 1 ],t; 



A + // - /i<9 tl log ■ 

77i(t + 2h[\- 1 } + 2%~ 1 ], t)77i(t,t) 
T*(i + 2h[X~ 1 ], t)r h (t + 2%^], t) ' 



(5.2) 



This equation contains three types of "shifted" r-functions. We can use the already mentioned 
differential operator 

oo _2 n _l 

D(z)=D(t,z)=J2^—7d t2n+1 

n=0 

to express these shifted r-functions as 

T h (t + 2hl\- l ],i) = e 2hD ^T h (t,i), 

m(t + 2h[^ 1 },i) = e 2nD ^m(t,t), 

r h {t + 2h[\~ l ] + 2% _1 ], t) = e 2hD ^ +2hD ^r h (t,i). 
This enables us to evaluate the logarithm of "r-quotients" in the previous equation as 

l0g ^2%-i]Jj = {e2hDW ~ e2HDM) l0gTa( *' * } = 2?l ~ 1(Z?(A) " Difl))F + ° (1) ' 
log ^(t + 2fe[A-|] + 2%-i],t) = (e2RD(A)+aftD&l) _ 1} logrn(t?i) 



log 



= 2/r 1 (D(A)+ J D(^)) J F + 0(l), 

r fi (t + 2^[A~ X ] + 2%~ 1 ], t)77i(t,t) 
T)i(t + 2ft[A- 1 ], i)r h (t + 2%^], t) 

= ( e 2KD(X)+2HD M + 1 _ e 2RD(A) _ ^M) logTfi ( t) t) = AD{\)D{n)F + 0{K). 



Having these results, we can readily find the limit of each terms in Q5.2[) as h — > 0. This yields 
the equation 

^(A -ix- 2D(X)d tl F + 2D^i)d tl F) = (A + n - 2D{X)d tl F - 2D(/j,)d tl F)e AD ^ D ^ F 
A — // 



for the F-function. We can rewrite this equation as 

(5.3) 



p(X) -p(fl) _ A-£ 4D(\)D(v.)F 



p(X)+p(fi) X + fj, 
where p(z) denotes the generating function 

°° 2z~ 2n ~ 1 
p(z) = z - 2D(z)d tl F = z-^2 -^pr dt2 " +AlF 

n=0 
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of dt 2n+1 dt 1 F. 1)5.3)1 is exactly the dispersionless Hirota equation derived by Bogdanov and 
Konopelchenko ^J for the 1-BKP hierarchy. We can also introduce the generating function 

oo 
S(z)=Y,t2n + lZ 2n+1 -D(z)F 
n=Q 

of the first derivatives of F to express p(z) as 

p(z) = d tl S(z) 

and rewrite (|5.3[) as 

d tl S(X) - d tl S([i) _ A - jx ad(X)d(u)f 
d tl S(\) + d tl S(v) \ + fi 

2) We can start from 1)4.2)1 in the rescaled form 
A + A*/\ h * i r h (t,t + 2h[X- 1 }) \ 

, Tft(t , t + 2ft[A~ 1 ] + 2%-i] )77i(t, t) 
X jfcCt, t + 2^[A- 1 ])r^(t, t + 2H\j,- 1 ]) 

and repeat the same calculations as above. The outcome is the dispersionless Hirota equation 



p(X)-p(fi) = X-£ e 4D{X)DQM)F 

of the 1-BKP hierarchy in the t-sector. D{z) denotes the differential operator 



(5.4) 



— 7 — in— i 

D(z) = D(i,z) =V 

V ' y ' ; ^ 2ra + 1 



' z~ 2n ~ x 

dt 2 n + l 



n=0 
and p(z) the generating function 

°° 2z~ 2n ~ 1 
p{z) = z - 2D{z)d h F = z-Y, -^T^dhF 

n=Q 

of c% 2? d^F. We can also use the generating function 

oo 

S{z) = Y J hn + iz 2n+1 - D{z)F 

n=0 

of the first derivatives of F to express p(z) as 

p(z) = d h S(z) 
and to rewrite 1)5 .4|) as 

%£(^) ~ %£Uf) _ A — // 4D(A)D(/i)F 

S fl 5(A)+a t - 1 5(/i) A + // 
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3) Let us now turn to (|4.3|) . In the rescaled form, it reads 

T h {t + 2h[\~ l lt) 



A - hd tl log 



r h (t, i + 2%- 1 ]) 

r^(t + 2h[\- 1 }, t + 2%-!]) \ r h (t + 2»[A- 1 ] J t + 2%- 1 ])r /i (t, t) 



A - fr<9 tl log ■ 



nit, t) J T h {t + 2h[\~ 1 }, t)m(t + 2%-i] , t) 



The logarithm of r-quotients in this equation can be evaluated in much the same way as the 
previous two cases. We thus obtain the equation 

A - 2D{X)d tl F + 2D(n)d tl F = (A - 2D(X)d tl F - 2D(fi)d tl F)e 4D ^ D ^ F 

for the F-function. In terms of the S-functions, this equation reads 

dhS(\) - d tl S(/J,) _ o 4d(\)d(jj,)F 



d tl S(X) + d tl S(n) 



e^WW*. (5.5) 



Unlike the previous two equations, this equation has no counterpart in the 1-BKP hierarchy. 
4) We can repeat the same calculations for (|4.4|) in the rescaled form 

T h (t, t + 2%- 1 ]) 
fi - hd h log r n - 

T h (t + 2h[\ i J, t) 

T h (t + 2h[\- 1 }, t + 2%- 1 ])^ T h (t + 2h[\~ 1 }, t + 2%- 1 ])r ft (t,t) 



// - Hdt 1 log 



r n (t,i) J r h (t + 2h[\~ 1 ], i)T h (t + 2%~i], t) 

to obtain the equation 

d^Sifl) - 5^ 5(A) _ o 4£>(A)D(/i)F 



%5( M )+%5(A) 



e ^Wiw_ (5. 6 ) 



This equation, like (|5.5|1 . has no counterpart in the 1-BKP hierarchy. 

As a final remark, let us note that the right hand sides of Q5.5JI and (|5.6|) coincide. This 
implies that the extra equation 

d tl S(X) - dtjjfj.) _ d h S(fi) - d h S(X) 
d tl S(X) + d tl S(fi) d k S(fi) + d h S(X) 

or, equivalently, 

d tl s(\) ■ d- tl s(x) = dtM ■ d h s(fi) 

holds. Since A and [i take arbitrary values, both sides of the last equation are actually indepen- 
dent of A and /i. This value can be determined by letting A, (X — > oo: 

lim d tl S(X) ■ d h S(X) = lim d tl S(n) ■ d h S(fi) = -2d tl d h F 

A-+oo £t— >oo 

Thus the extra equation reduces to 

d tl S(z) ■ d h S(z) = d tl S(z) ■ d h S(z) = u, (5.7) 

where 

u = -2d tl d h F. 

As it will turn out later, this equation is the Hamilton-Jacobi equation of a two-dimensional 
Schrodinger equation that underlies the Novikov-Veselov hierarchy [23J. 
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6 Hamilton— Jacobi equations 

We show here that the dispersionless Hirota equations (|5.3|) - (|5.6|) are equivalent to a system of 
Hamilton-Jacobi equations for the ^-functions. As we shall show later, these Hamilton-Jacobi 
equations are Lax equations of the dispersionless 2-BKP hierarchy in disguise. In other words, 
this system, too, is a special case of the universal Whitham hierarchy |16j . 

Following Teo's idea [2^, we borrow the notion of "Faber polynomials" from the theory of 
univalent functions (see Teo's paper [2J and references cited therein). This enables us to do 
calculations in a clear and systematic way. The Faber polynomials are defined by a generating 
function, and eventually turn out to give the Hamiltonians of the Hamilton-Jacobi equations. 

1) To derive Hamilton-Jacobi equations from (J5.3J) . we first rewrite it to the logarithmic 
form 

We now introduce the Faber polynomials $ n (p), n = 1,2,..., of p(z). They are defined by 
a generating function of the form 

viz] — w v-^ z~ n T . . . 

log — = -V $ n (w). (6.2) 

z ^-^ n 

n=l 

Although this definition is valid for any Laurent series p(z) of the form p{z) = z + 0(1), the 
present setting is special in that p(z) is an odd function: 

p(-z) = -p{z). 

Because of this, we have the identity 

p(— z)+w , p(z)—w 
log = log , 

— z z 

which implies that & n (w) is even if n is even and odd if n is odd, i.e., 

$n(-w) = (-l) n $nM- 

Therefore we have 

log E^—^ = -J2 — ($„H - *„(-«,)) = - £ ?f^<S> 2n+1 ( w ). (6.3) 

v[z)+w t—* n ^-^ 2n + 1 

^ v ; n=l n=0 



Substituting z = A and w = p{a) yields 

PW+PW ^ 2n+l 

which coincides with the left hand side of (|6.1|) . As regards the right hand side of (|6.1|) . we can 
rewrite the double sum as 

00 4^-2n-l -2m-l °° 2\~ 2n ~ 1 ( °° 2u~ 2m_1 \ 

^ (2n + l)(2m + l) dt2n+ldt2m+lF = 5-, 2n + 1 ( 2 2 m + 1 9 W 9 Wh^J 

m,n=0 y Jy ' n=0 \m=0 / 

00 o\-2n-l \ _ ,, °° o\-2n-l 

= y t± u2« + i _ d s{ }) = _ Y d_j±-y?£ 9 S (n). 

*-" 2n+l VP t2n+1 VPy; A + u Z^ 2n + 1 * 2n+1 VPy 

n=0 n=0 
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Thus 1)6.1)1 reduces to the following system of Hamilton-Jacobi equations: 

d t2n+1 S(z)=$2n+l(p(z)) = $2n+l(dt 1 S(z)), 71 = 0,1,.... (6.4) 

2) Let us now consider Q5.4)) . In this case, we use the Farber polynomials $ n (p), n = 1,2,... 
of p(z) defined by the generating functional relation 

. p(z) -w ^ z~ n - 

log ^-^ = - \ $ n (w . 

z *-^ n 

n=l 

Since p(z), like p(z), is an odd function of z, we have the generating function 

. p(z)-w ^ 2z ~ 2n ~ 1 ^ i S 

pm + it; ^— ' 2n + 1 

v ' n=0 

of the Faber polynomials with odd indices. The rest of calculations are almost the same as in 
the case of 1). Thus 1)5. 4JI reduces to the following Hamilton-Jacobi equations: 

dt 2n+1 S(z) = $2n+i(p(z)) = $2n+i{d h S{z)), n = 0, 1, ... . (6.5) 

3) We can deal with Q5.5)) in much the same way as the case of 1). We start from the 
logarithmic form 

p(A)-Mfr) y 4A- 2 -V 2 -i 
g p(A) + MM ^ (2n + l)(2m + 1) * 2 " +1 t2m+1 

of (|5.5|) . By 1)6.3)1 . the left hand side can be expressed as 



The right hand side can be written as 

00 A\—2n—l..—2m—l ^ n\- 

5- (2n+l)(2m+l) a42 " +iaf2m + lF = " 2 "^ f^Wi* W 

n,m=0 n=0 

Thus 1)5.5)1 reduces to Hamilton-Jacobi equations of the form 

9t 2n+1 S(z) = $2n+i(d tl S(z)), n = 0,l,.... (6.6) 

4) ()5.6|) can be treated in the same way as (|5.4[) . and boils down to the following Hamilton- 
Jacobi equations: 

a f2n+1 S(z)=$ 2n+1 (a 4 - 1 5(z)), n = 0,l,2,.... (6.7) 

We have thus derived the Hamilton-Jacobi equations 1)6. 4J) . 1)6. 5)1 . 1)6.6)1 and 1)6.7)1 from the 
dispersionless Hirota equations 1)5.3)1 . ()5.4)l . 1)5.5)1 and 1)5.6)1 . Since what we have done is simply 
to expand generating functions, this procedure can be reversed. The equivalence of these two 
systems are thus proven. 

As a final remark, we note that these four sets of Hamilton-Jacobi equations can be packed 
into total differential equations of the form 

oo oo 

dS{z) = S'(z)dz + Y,$2n + i(d tl S(z))dt 2n+1 +J2^2n + i(dt 1 S(z))di 2 n + i, (6.8) 

n=0 ra=0 
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oo oo 

dS{z) = S'{z)dz + Y,$2n + i(d tl S(z))dt 2n+1 +Y / $2n + i(di 1 S(z))di 2n+l , (6.9) 

n=0 n=0 

where S'(z) and S'(z) denote the z-derivative of S(z) and S(z), 
S'(z) = ^1, S'(z) - d ^ z) 



dz dz 

7 Dispersionless Lax equations 

The Hamilton-Jacobi equations (|6.4|) . (|6.5|) . (|6.6|) and (|6.7j) yield a system of "dispersionless Lax 
equations" (i.e., Lax equations with respect to Poisson brackets rather than usual commutators) 
for the inverse functions C(p), C{p) of p(z), p(z). Since the derivation is almost the same as 
the case of the dispersionless KP and Toda hierarchies [TJI], we here show just an outline of the 
results, specifying some new aspects. 

7.1 How to derive dispersionless Lax equations 

Let z = C(p) denote the inverse function of p = p(z) in a neighborhood of z = oo. This is an 
odd function and has Laurent expansion of the form 

dp) =p+o{p- 1 ) 

in a neighborhood of p = oo. We now substitute z = C(p) into (|6.8|) and note that ()5.7j) implies 
the relation 



a k s( z y 



=Ap) d tl S(z) 
This yields the equation 



u 

z=C(p) P 



dS{p) =M{p)dC{p) +pdh + ^2^2n+l(p)dt 2 n + l + -&1 + J^ W(")4^1, (7-1) 

n=l n=l 

where 

S(p) = S(C(p)), M{p) = S'{C{p)). 

Similarly, we can consider the inverse function z = C{p) of p = p(z), which is an odd function 
of p and has Laurent expansion of the form 

C(p) =p + o{p- 1 ) 

in a neighborhood of p = oo. Since Q5.7JI implies that 



d tl S(z)' 



-m d h s(z) 

H6.9J) turns into the equation 



z=C(p) 



u 

P 



dS(p) =M(p)dC(p) + -dtl + J2®2n+l(-)dt2n+l + pdtl + ^ $2n+l(p)d*2n+l, (7.2) 



n=l n=l 
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where 

S(p) = S(£(p)), M(p) = S'(C(p)). 

C(p), L{p) and A4(p), M(p) are dispersionless analogues of the Lax and Orlov-Schulman ope- 
rators of the KP and Toda hierarchies jlUj . Actually, as we show below, the present situation 
is more complicated than the case of the KP and Toda hierarchies because of presence of two 
"spatial variables" t\ and t\. 

p and p in 1)7.1(1 and (|7.2j) are a priori independent. To derive dispersionless Lax equations, 
we have to set the algebraic relation 

pp = u 

between p and p. Note here that the Faber polynomials in the coefficients of (|7.1[) and ([7.2)1 
thereby coincide as 

$2n+l (p) = $2n+l (u/p) , <§2n+l («/*») = $2n+l (p) • 

Therefore we have the well-defined 1-form 

oo oo 

9 = pdt x + ^ ®2n+l(p)dt2n+l + pdt\ + ^J $2n+l(p) d *2n+l 
n=l n=l 

and the equation 

dS(p) - M(p)dC(p) = = dS(p) -M{p)dC{p) 
of 1-forms, from which we obtain the equation 

dC{p) A dM{p) = dC(p) A dM(p) 

of 2-forms. 

Having these equations, we can now resort to the standard procedure ^U] to derive disper- 
sionless Lax equations. Actually, because of the manifest symmetry between t, p, £(p), M.(p), 
<3?2n+i(p) and i, p, £(p), A-^(p), &2n+i(p), we can formulate the dispersionless Lax equations in 
two different ways. A formulation is based on the Poisson bracket 

if \ = ^£^L d f dff 
U>9h,ti dpdti dhdp 

for the canonical pair (p, ti). We have the dispersionless Lax equations 

dt 2n+ Mv) = {$2n+i(p),£(p)W, d hn+1 K{p) = {$ 2n+1 (u/p),lC( P )} pM (7.3) 

for IC(p) = C(p), A4(p), C(u/p), M.(u/p) and the canonical commutation relations 

{C(p),M( P )} p , tl = 1, {C(u/p),M(u/p)} p , tl = 1. 
Another formulation is based on the Poisson bracket 

If \ = ^L^L df dg 

U'dtpM dpdli dlidp - 

for the canonical pair (p, ii). We have the dispersionless Lax equations 

dt 2n+ MP) = {$2n + i{u/p),1C(p)}p,h, dt 2n+ MP) = {^n+imXmpM (7.4) 

for IC(p) = £(p), A4(p), £(u/p), M.(u/p) and the canonical commutation relations 
{£(p),M(p)}p fl = 1, {C(u/p),M(u/p)} pA = 1. 
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7.2 How to determine Hamiltonians 

For these dispersionless Lax equations to form a closed system, the Hamiltonians have to be 
written in terms of the L- functions. 

We can read off such an expression of the Hamiltonians from the Hamilton-Jacobi equations 
themselves. Let us examine (|6.4JI . Substituting z = C{p) into this equation yields 

<&2n+i(p) = d t2n+1 S(z)\ z=c{p) = C(p) 2n+1 + OiCip)- 1 ), 
which implies that <3?2n+i(p) coincides with the "polynomial part" of C(p) 2n+l : 

<$> 2n+l (p) = {C{p? n+l )^. (7.5) 

In the same way, (|6.7|) yields 

§ 2n+1 (p) = d hn+1 S(z)\ z=m = C(p) 2n+1 + Oip- 1 ), 
so that <&2n+i(p) is given by the polynomial part of the powers of C(p): 

<W(p) = (£(p) 2n+1 )> - (7-6) 

Actually, (|7.5|) and H7.6|) can be derived from the definition of the Faber polynomials as well. 
As regards <&2n+i(p)j the defining relation l|6.2|) differentiated by z yields the identity 



p'(z) 



xy^^n- 



P(z) ~ w n=1 

Therefore, by Cauchy's integral theorem, we have 
dz z n p'(z) I dp C(p) n 



$n(w) 



C 2m p(z) - w J\p\=R 2ttz p-w 



where C is the image of the circle \p\ = R by the map z = £(p), being a simple curve for a 
sufficiently large value of R, and w is understood to sit on the near side, i.e., |to| < R. If we 
insert the Laurent expansion of C(p) in this expression, we obtain the algebraic expression (|7.5|) . 
In particular, (J7.5I) holds for the Faber polynomials with even indices as well, though we do 
not use them here. Similarly, we have an integral formula of $ n (p), which implies the algebraic 
expression (|7.6|) and its counterpart for even indices. 

8 Another approach to dispersionless limit 

In this section, we present another approach that was taken in Takasaki's previous work on 
the 1-BKP hierarchy [Hj. This approach is based on quasi-classical (or "WKB") approximation 
of auxiliary linear equations (see the review JO] and references cited therein). As expected, 
this approach leads to the same Hamilton-Jacobi equations as derived from the differential Fay 
identities. 

8.1 Wave function and auxiliary linear equations 

Since we have two vertex operators X(t, z) and X(t, z), we can define two wave functions: 

T{t,t) T(t,t) 
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More explicitly, 



r(t - 2[z-\ t) Mt , z) = {f , , _ T(t,t-2[^]) ^ zl 



The bilinear equation (|2.1I) for the tau function thereby turns into the bilinear equation 

Z -#(*', t', *)tt(t, t, -*) = / ^r-Vtf, t', z)#(t, t, -z) (8.1) 



27riz / 27riz 

for the wave functions. Note that the role of "dual wave functions" in the KP and Toda 
hierarchies are now played by *$>(t, i, —z) and ^>(t, i, —z). Let us mention that such a pair of wave 
functions for the 2-BKP hierarchy were first considered by Shiota [7] in an algebro-geometric 
framework. 

As we shall show below, these two wave functions satisfy the following three sets of "auxiliary 
linear equations." 

1. Linear equations of the 1-BKP type in the t-sector: 

{d t2n+1 -B 2n+1 (d tl ))y = 0, n = 0,1,2,..., (8.2) 

where B2 n +i(dt 1 ) is a differential operator in t\ of the form d t n+ + • • • without 0-th order 
term, i.e., 

B 2n+1 (d tl )l = 0. 

2. Linear equations of the 1-BKP type in the t-sector: 

(% n+1 - B 2n+l (d h ))V = 0, n = 0, 1, 2, . . . , (8.3) 

where B 2n +i (d^ ) is a differential operator in t\ of the form d^ n+ + • • • without 0-th order 
term, i.e., 

B 2n+1 (d- h )l = 0. 

3. Two-dimensional Schrodinger equation of the Novikov-Veselov type |24j : 

(d tl d h -u)=0 (8.4) 

where u = — 2c?£ 1 c% 1 logr. 

Note that the condition on the 0-th order term of B 2n +i(dt 1 ) and B 2n +i(dt ,) is the same as 
the condition that characterizes the .B-operators of the 1-BKP hierarchy [1JH)]. This is another 
way to see that the 2-BKP hierarchy contains two copies of the 1-BKP hierarchy. 

8.2 How to derive auxiliary linear equations 

These auxiliary linear equations are derived from the bilinear equation (|8.1|) by the standard 
method H HI H3 |U] based on the following 

Lemma 1. If R(t, i, z) and R(t, i, z) are a pair of Laurent series of z of the form 
R(t, t, z) = 0(z- V ( H R(t, t, z) = 0(z" V ( H 

and satisfy the bilinear equation 

f ^k R{t '' *'' z)q,{t > *' ~ z) = f ^k R{t '' *'' z) * (t> *' ~ z) ' 

then R(t, i, z) = R(t, i, z) = 0. 
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We refer the proof of this lemma to the references 1221 1291 I3UJ and illustrate usage of this 
lemma for the case of (|8.2|) . Differentiating (|8.1|) by t 2n+ \ yields 

/dz - I" dz 
d t i ^(t', *', z)W(t, t, -z) = & A, *(t' *', z)tf (t, i, -A 

Similar bilinear equations hold for higher derivatives as well. In particular, we have 

j> ^B(d tl )V(t', t', z)*(t, t, -z) = <f> ^L B (d tl )*(t', t', z)*{t, t, -z). 

for any choice of differential operator B{dt 1 ) in t\. This implies that the bilinear equation in 
the lemma is satisfied by 

R(t, i, z) = d w *(£, i, z) - B(d tl )V(t, i, z) , 
R(t,i,z)=d t2n+I y(t,i,z)-B(d tl )y(t,i,z). 

If B{dt 1 ) is chosen to be an operator B2 n +i{dt 1 ) for which other conditions of the lemma are 
satisfied, then R(t,i,z) and R(t,t,z) automatically vanish, thus we obtain (|8.2j) . The extra 
condition on the 0-th order term of B2 n +i(dt 1 ) turns out to be satisfied because (|8.1|) reduces 
to the bilinear equation of the 1-BKP hierarchy in the t-sector upon setting i =i. 
(|8.3|) and (|8.4j) can be derived in the same way. 

8.3 Lax and zero-curvature equations 

We now consider a "scalar" Lax formalism of the 2-BKP hierarchy on the basis of the auxiliary 
linear equations presented above. This Lax formalism is conceptually parallel to the algebro- 
geometric framework based on "two-point" (or "two-puncture") Baker- Akhiezer functions [7[|Sj. 
In spite of some unusual aspects as we see below, such a scalar Lax formalism is more suited 
to the consideration of dispersionless limit than the "matrix" Lax formalism of Kac and van de 
Leur [22] . 

Let us first consider zero-curvature equations for the .B-operators. As pointed out by Kri- 
chever |5], the S-operators satisfy the modified "zero-curvature" equations 

[dt 2m+1 - B 2m +i(dt 1 ), d hn+1 - B2 n+ i{d h )] = D mn (d t ,di t ){d tl di t - u), 

where D mn (dt 1 ,di 1 ) are differential operators in both t\ and i±, alongside the usual zero- 
curvature equations 

[ d t 2m +i ~ B 2m +i{d tl ), d t2n+1 - B 2n +i{d tl )] = 0, 
[ d hm+i ~ B 2m +i(d- tl ), d hn+1 - B 2n+1 (d h )] = 

of the 1-BKP hierarchy in the t- and t-sectors. In other words, c\ 2m+1 — B 2m+ \{dt 1 ) and dj 1 — 
B 2n +i(dt 1 ) commute "modulo" the two-dimensional Schrodinger operator dt 1 dt 1 — u. We refer 
the details to Krichever's paper [H]. 

Let us turn to Lax equations. We have two L-operators 

L(d tl ) = w(d tl )d tl w(d tl r\ L(d h ) = w(d h )d h w(d- tl r\ 

where W^d^) and Wi^d^) are dressing operators of the form 

oo oo 

W(d tl ) = 1 + Y,™j(t,i)d^, W(d h ) = l + ^^(t,t)^' 

3=1 3=1 
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that are related to the wave functions as 

V(t,i,z) = W(d tl )e^ z \ *(t,t,z) = W{d h )e^ l ' z \ 

Substituting this expression of the wave functions into (|8.2j) and (|8.3j) yields the so called Sato 
equations 

^^ = B 2n+l (d tl )W(d tl ) - W{d tl )d^\ 

dw{d h 



B 2n+l (d h )W{d h )-W(d h )df + \ 



9t2n+l 
which in turn imply the following expression of the 5-operators: 

B 2n+l (d tl ) = (L(d tl ) 2n+1 )> Q , B 2n+l {d~ tl ) = (L(d h ) 2n+1 )> , (8.5) 

where ( )>o denotes the differential operator part of a pseudo-differential operator. Since 
^(t,i,z) and ^f(t,i,z) satisfy the linear equations 

L(d tl )W(t,i,z)=zW(t,i,z), L(t,i,z)V(t,i,z) = z#(t,t,z) (8.6) 

alongside l|8.2|l and (|8.3j) . we have the Lax equations 

^^ = [B 2n+1 (d tl ),L(d tl )], ^^ = [B 2n+1 (d h ),Hd h )], (8.7) 

which give the standard Lax representation of the 1-BKP hierarchy in the t- and t-sector, 
respectively. 

It is, however, not so straightforward to derive Lax equations of the form 

dL(d tl ) dl{d h ) - 

-= = [*,L(d tl )\, — — '- = [*, L(d h )}, 

ot 2n+ i ot 2n +i 

Those Lax equations can be obtained from auxiliary linear equations of the form 

{Q2n + l(d tl )dt 2n+1 - P 2n+1 (d tl ))V = 0, {Q2n+i(dt 1 )d t2n+1 - P 2n+1 (d h ))^> = 0, (8.8) 

where P 2n+ i(d tl ), Q 2n +i(9t 1 ) and P 2n +i(dt 1 ), Q 2n +i(h) are differential operators in t\ and ii, 
respectively. One can derive these linear equations by the same reasoning as in the case of (|8.4(l . 
Note that (|8.4|) amounts to the case of n = 1 with 

Qi (d tl ) = d tl , Pi (d tl ) = u, Qt (% ) = % , Pi (d h ) = u. 

Since (|8.6j) and (|8.8() are compatible, the L-operators turn out to satisfy the Lax equations 

^£^ = [Q2n + l{d tl Y l P2n + l{d tl ), !(%)], 



Of 



2n+l 



^*± = [Q 2n+1 (d h )- 1 P 2n+1 (d h ), L(d h )]. (8.9) 



Dispersionless Hirota Equations of Two-Component BKP Hierarchy 19 

8.4 Quasi-classical approximation of wave functions 

We now consider the limit as h — » 0. If the rescaled r-function Tft(t,i) behaves like (|5.1|) . the 
associated wave functions 

turn out to take the quasi-classical form 

* R (t, t, z) = exp(fT l S(t, t, z) + 0(1)) , f fi (t, t, z) = exp(/T ^(i, i, z) + O(l)) . (8.10) 

To simplify notations, let us write the phase functions S(t,i,z), S(t,i,z) as S(z),S(z). As 
expected, they have the Laurent expansion 

00 °° o„-2n-l 

n=0 n=0 

00 °° o„-2n-l 

s(z) = E w^ +1 - £ |r+r^^ 

n=0 n=0 

in a neighborhood of z = oo. As the time variables are rescaled, differential operators in the 
auxiliary linear equations are also rescaled as 

Substituting the quasi-classical form (|8.1U|) of the wave functions into those linear equations, we 
obtain a set of Hamilton-Jacobi equations for the phase functions. 

The fundamental auxiliary linear equations (|8.2|) . (|8.3j) and (|8.4|) thus yield the Hamilton- 
Jacobi equations 

dt 2n+ iS{z) = B 2n+1 (d tl S(z)), 
d hn+1 S(z)=B 2n+1 {d h S(z)), 
d tl S(z)-d h S(z) = u (8.11) 

for S(z) and 

d t2n+1 S(z)=B 2n+1 {d tl S(z)), 
d hn+1 S(z)=B 2n+1 {d h S(z)), 
d tl S(z) ■ d h S(z) = u 

for S(z). B(p) and B(p) denote the classical limit 

#2n+i(i>) = hm B 2n +i(p), B 2n+ i(p) = lim B 2n+1 (p) 
h-^o h— >o 

of the "symbols" of the corresponding differential operators. 
Similarly 1)8.6(1 yields the equations 

C(d tl S(z)) = z, C(d h S(z))=z, 

where C(p) and C(p), like B 2n+ i(p) and B 2n+ i(p), are given by 

C(p) = limL(p), C(p) = lira L(p). 

h^o fi,— >o 
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Thus p(z) = dt 1 S(z) and p(z) = dt\S(z) turn out to be the inverse functions of z = C(jp) and 
z = C(p), respectively. The relation 1)8.5)) connecting the S-operators and the L-operators turns 
into the relation 

B 2n+l (p) = (£(p) 2 " +1 )> , B 2n+1 (p) = (C(p) 2n+1 )> 

among the corresponding phase space functions. (-)>0 now stand for the polynomial part of 
Laurent series (of p and p, respectively). This, in particular, implies that B 2n+ i(p) coincide with 
the Farber polynomials of the inverse function p = p(z) of z = C(p). A similar statement holds 
for B 2n +i{p) and the inverse function p = p(z) of z = C{p). Thus we have 

B 2n +l{p) = $2n+l(p), B~2n+l{p) = $2n+l(p)- 

The somewhat strange linear equations (|8.8|) . too, turn out to have a natural interpretation. 
The associated Hamilton-Jacobi equations read 

Q 2 n+i(d tl S(z))d hn+1 S(z) = p2n+i(d h S(z)), 

Q 2n+1 (d h S(z))d t2n+1 S(z) = V 2n+1 (d h S(z)), 

for S(z) and 

Q 2n+l (d tl S(z))d hn+1 S(z) = P 2n+l (d tl S(z)), 
Q 2n+l {d- tl S(z))d t2n+1 S(z) = V 2n+1 (d h S(z)), 

S(z), where 

Q 2n +i{p) = lim Q 2n+ i(p), P 2n+ i(p) = lim P 2 n+i(p), 

Q.2n+i(P) = lim Q 2n+ i(p), T> 2n +i(p) = lim P 2n+ i{p). 
h^o h— >o 

For the equations of the t-flows to be consistent with the previous Hamilton-Jacobi equations, 
V 2n+ i(p), Q 2n +i{p) and B 2n+ \{p) have to satisfy the equations 

Q 2n+ i(d tl s(z)y 1 p 2n+1 (d tl s(z)) = B 2n+1 (d h s(z)), 

Q2n + i{d tl S{z)Y l P 2n+l {d tl S{z)) = B 2n+l {d h S(z)). 

Since dt 1 S(z) and di 1 S{z) obey the algebraic relation (|8.11|) . these two equations reduce to a 
single equation of the form 

B 2n+1 {u/p) = Q2n+l{py 1 p2n+l(p)- (8.12) 

Similarly, we have the relation 

B 2n+1 (u/p) = Q 2 n+l(p)" 1 ^2n+l(p). (8.13) 

Actually, this is a result to be expected from 1)8. 8JI . because 1)8. 8jl can be (formally) rewritten as 

{d hn+1 ~ Q2n + l(dt i r 1 P2n+l(dt 1 ))* = 0, 
[d t2n+1 ~ Q2n + l(d h )- l P2n+l(dt 1 ))* = 0, 

and for these equations to be consistent with ()8.2|) and 1)8. 3|) . we should have 

B 2n+l (d- tl )^> = Q 2n+l {d tl )- l P 2n+l (d tl )^, 

B 2n+1 (d tl )V = Q 2n+l (d h )- 1 P 2n+1 (d h )y, 
which may be thought of as "quantization" of 1)8.12)1 and (|8.13)l . 
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8.5 Dispersionless Lax equations 

In the limit as H — > 0, (|8.7|) and (|8.9|) turn into the dispersionless Lax equations 

d t2n+i£(p) = {#2n+l(p), C{p)} pM , d hn+1 C(p) = {B 2 n+l(p)X(P)}p,t 1 

and 

dt 2n+1 C(p) = {Q 2 n+i(p)" 1 i 3 2n+i(p), C{p)} pM , 

d t2n+l£(P) = {Q2n+l(p)~ 1 'P2n+l(p), C{p)} pfl , 

respectively. In view of 1)8.12(1 and (|8.13|) . we can rewrite the second two as 

<9f 2n+1 £0) = {B2n+i(u/p),£(p)} p , tl , d t2n+1 £(p) = {B 2 n+i(u/p),C(P)}p,t 1 - 
Obviously, these equations coincide with the corresponding equations in (|7.3|) and (|7.4|) , 
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